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Abstract. The peculiar motion of an observer with respect to the CMB rest frame induces 
an apparent deflection of the observed CMB photons, i.e. aberration, and a shift in their 
frequency, i.e. Doppler effect. Both effects distort the temperature multipoles a£ m 's via 
a mixing matrix at any I. The common lore when performing a CMB based cosmological 
parameter estimation is to consider that Doppler affects only the I = 1 multipole, and 
neglect any other corrections. In this paper we reconsider the validity of this assumption, 
showing that it is actually not robust when sky cuts are included to model CMB foreground 
contaminations. Assuming a simple fiducial cosmological model with five parameters, we 
simulated CMB temperature maps of the sky in a WMAP-like and in a Planck-like experiment 
and added aberration and Doppler effects to the maps. We then analyzed with a MCMC in 
a Bayesian framework the maps with and without aberration and Doppler effects in order 
to assess the ability of reconstructing the parameters of the fiducial model. We find that, 
depending on the specific realization of the simulated data, the parameters can be biased 
up to one standard deviation for WMAP and almost two standard deviations for Planck. 
Therefore we conclude that in general it is not a solid assumption to neglect aberration in a 
CMB based cosmological parameter estimation. 
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1 Introduction 

The Cosmic Microwave Background is used as a fundamental tool to test Cosmological Mod- 
els and to quantitatively extract the parameters of such models. This is usually done by 
extracting the Temperature and the Polarization power spectra from the maps and by fitting 
them with a ACDM background model with an almost scale-invariant gaussian spectrum of 
density perturbations. However, if we observe the CMB with a large velocity /3 = v/c relative 
to such background, the image undergoes distortions due to the Doppler effect and due to 
the aberration effect and it is natural to wonder whether there is any sizable bias on such 
spectra and so also on cosmological parameter estimation due to the fact that we are not at 
rest with the CMB. 

It is always commonly assumed that the only sizable effect of an observer's velocity on 
the CMB map is the generation of a large dipole and any other possible effect is commonly 
ignored in all analyses. In fact, if we consider a perfectly homogeneous map, there is a 
Doppler effect which induces a dipole of order (3, while higher n-th multipoles are of order f3 n . 
Since we observe a dipole of order 10 -3 , while other multipoles are 0(1(T 5 ), the standard 
lore is to interpret this large dipole as due to our velocity (which is then inferred to be 
P = (1.23 ± 0.003) x 10~ 3 [1]), and otherwise to ignore completely any other effect of our 
velocity. 

However this would be true only on a perfectly homogenous sky. In reality we have a 
nontrivial map with nonzero coefficients ag m , which are measured to be of order 10 -5 for any 
i > 1. Applying a Doppler effect induces a mixing between multipoles [2-4]. Precisely there 
is a mixing at order /3 between neighboring multipoles I and I db 1 and a mixing at order j3 n 
among I and £±n. Moreover there is another important distortion, that is aberration, which 
changes the direction of observation of incoming photons, similarly to what lensing does. 
Aberration [2-4] can be computed for small j3 and induces an additional coupling between 
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neighboring multipoles of order (/3£) n . Such a distortion is therefore even more important 
because it becomes large at high I. Note that this is in fact the most important secondary 
effect on the CMB and it is not clear a priori if it induces a negligible bias on the power 
spectra and on the cosmological parameters: an 0((3) or 0(f3£) effect on all multipoles could 
a priori lead to a large effect, when summed over all £ and rrt's. Moreover for aberration this 
effect grows as (3£ and becomes dominant at large t: using the value (3 = 1.23 X 1CT 3 the 
distortion is order 1 on a single coefficient already at £ > 800! What happens in fact is that 
the aberration is deflecting the arrival direction of the photons by an angle of order (3 and so 
if we are looking at angles smaller than f3 there is an order 1 effect on the CMB. 

However, a negative result came from [2] where it was found that, defining as usual C% 
I [ 2 I + 1) Ylm a h ffl <mi t ne bias on the ensemble average (Ci) undergoes some cancellations 
leading to a very small deviation which does not grow with £, precisely 5(Ci) / (Ci) = 0(/3 2 ) 
at leading order in (3. 

More recently instead it has been realized that such conclusion is not the final word for 
several reasons. First of all such calculations are based on a perturbative approach in (3 and 
so, as we mentioned before, they apply only up to £ < 800, while we do not know if a clear bias 
could show up at higher ts [4, 5]. Second, such an estimate is based on a full-sky definition 
for Ci : in fact, in full-sky, because of the sum over m there are cancellations between the 
corrections on the individual ai m 7 s. However we never observe a full-sky in an experiment 
and it can be easily seen that as long as one part of the sky is removed such cancellations 
do not hold and we can get again effects of order f3£, which are larger, on the pseudo-Q's 
which have to be used in a partial sky estimate. This was noticed in [6], where the bias on 
the pseudo-C^s has been analyzed for small j3£ and for very low £, less than about 20. This 
estimate becomes very expensive computationally at higher £ and anyway not feasible within 
pertubation theory at large £ > 1/(3 . For all these reasons the question whether the effect 
of velocity is negligible on the power spectrum is still lacking an answer. 

The scope of the present paper is to resolve this issue by analyzing the bias on cosmo- 
logical parameter estimation due to a nonzero (3 at any £ and for both full-sky and cut-sky. 
As we have stressed this could be done perturbatively analytically only for \ow-£, while an- 
other possibility, valid at any £, would be to write down the exact mixing matrix between 
any multipole without pertubative expansions. This would amount to numerically compute 
integrals of highly oscillating Legendre polynomials, which becomes very slow at high-^, and 
becomes extremely slow for a cut-sky, since it involves also a window function. This approach 
has been addressed by recursive formulae in [7] and by approximate fitting functions in [5]. 

The way we proceed here is a more direct and straightforward one: we simulate maps of 
the CMB sky and we directly apply on the maps the boost trasformation before extracting 
the afrn's, therefore bypassing the need of computing the mixing coefficients. Then we extract 
the aim's and the C^s and finally we run a Markov Chain Monte Carlo (MCMC) to estimate 
cosmological parameters and evaluate the bias between the moving observer and the rest 
frame observer. This approach also allows us to easily implement sky-cuts and to check 
whether the effect is larger on a realistic cut-sky experiment. We apply all this procedure to 
Temperature maps only and both for WMAP-like resolution and for Planck-like resolution. 

The paper is organized as follows: in section 2 we review the impact of aberration 
and Doppler on CMB observables. Section 3 is instead devoted to a detailed description of 
the procedure followed to simulate CMB temperature maps properly including Doppler and 
aberration effects. These maps are then analyzed in section 4, while section 5 summarizes 
our results. 
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2 CMB aberration and Doppler 



For an observer in motion with respect to a light source, the aberration phenomenon consists 
in the apparent deflection of the observed light bundles due to the relative motion between 
the observer and the source. In the case of the cosmic microwave background (CMB), as a 
consequence of aberration, an observer moving with respect to the CMB rest frame assigns 
to a photon emitted along the direction n at the last scattering surface the arrival direction 
n', where n and n' are unit vectors. The aberration angle a = h' — h can be calculated from 
the velocity transformation relating the CMB and the observer reference frames: assuming 
that the observer moves along the direction of the z-axis (identified by the unit vector z) 
with speed /3 compared to the CMB reference frame, one finds 

A P sin 2 6 

a-z = (2.1) 

l + /3cos# V ; 

where cos# = n • z. We refer to the relative motion between the observer and the CMB rest 
frame as the peculiar motion of the observer. Besides an aberration of the CMB radiation, 
a peculiar motion of the observer induces also a change in the frequency of the observed 
photons. According to Special Relativity the frequency v' in the moving frame is related to 
the frequency is measured at rest by the usual Doppler effect: 

v' = 1/7(1 + /3n ■ z) (2.2) 

where 7 = 1/ yl — /3 2 . Doppler and aberration effects unavoidably distort the CMB tem- 
perature and polarization maps currently measured by experiments such as for instance the 
WMAP or Planck satellites. This can be seen as follows: CMB experiments directly measure 
the brightness of the observed radiation. The brightness I (is) is then conventionally trans- 
lated into an equivalent value for the thermodynamic temperature T through the black body 
law 

I(is) e v ' T - 1 
IS is/ 1 

(2.4) 



By solving Eq. (2.3) one finds 



log(l+ ' 



which is therefore the temperature indirectly measured by CMB experiments. Analogously, 
in the case of the CMB polarization one translates the quantities actually measured into 
the value of an effective temperature obtained replacing in Eq. (2.4) the brightness I (is) 
with the relevant Stokes parameters. Since the ratio is 3 /I(is) is invariant with respect to 
frame transformations [2], one can conclude from Eq. (2.4) that under a Lorentz boost the 
temperature has to transform as follows 

T'(n') = 7(1 + /3n • z)T(n) (2.5) 

where T" is the temperature in the moving frame and T is the temperature measured by an 
observer at rest (with respect to the CMB). Eq. (2.5) together with Eq. (2.1) thus describes 
how a temperature map is distorted consequently to aberration and Doppler effects: the 
temperature associated with the direction n' by an observer in motion with peculiar velocity 
/3z is equal to the temperature that an observer at rest would assign to the direction n, mul- 
tiplied by a corrective Doppler factor which depends on (3. Expanding in spherical harmonics 
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both T'(n') = Ylim a 'im Y lm(n') and T(n) = Y l im a lmXlm{^), Eq. (2.5) can be rewritten as 
follows 1 

a i'm' = diia em> [t(1 + /5n • z)p 2 y^ m ,(n)l£m'(n) • ( 2 - 6 ) 

This expression can be then used to compute the desired correlation functions and com- 
pare these results with observations, properly accounting for aberration and Doppler effects. 
Eq. (2.6) is exact, though it relies on an integral computationally very expensive for high 
multipoles. 

In the following subsections we will review two approximations recently discussed in 
the literature to estimate the coefficients a' lm of Eq. (2.6). In sections 3 and 4, instead, 
we will show how starting from Eq. (2.5) one can actually include aberration and Doppler 
effects exactly in a CMB based cosmological parameter estimation, by transforming directly 
simulated maps of the sky in the boosted frame. 

2.1 Perturbative expansions 

Since our peculiar velocity is small in natural units, i.e. (3 ~ 1.23 x 1CP 3 , Eq. (2.6) can 
be evaluated expanding the integrand around (3 equal to zero. The leading corrections in (3 
necessary to evaluate the C/s are given by 



a 'im — c em ae-i m + cj m a e+lm + a im (l + d im ) (2.7) 



where 



-WGt 
2 



[(£ + i)(e + 2)Gj +hm + e(e-i)Gl m -e(£ + i) + m 2 -i] (2. 



and Gi^ m = J • Similar expressions can be derived for the analogous coefficients cor- 
responding to the Stokes parameters describing the CMB polarization. We refer the reader 
to Refs. [2-5] for an explicit derivation of such expressions; note that the correct expressions 
for Temperature fluctuations for the cf m have been derived first in [3] and for d£ m in [5] and 
they slightly differ from [2], which would be valid for Intensity instead of Temperature. We 
note here that the corrections increase as powers of f3£, and so become large for high £ and 
so for £ > 1/(3 the perturbative approach breaks down. This happens because the aberration 
effect is a distortion on the angle of order 1/(3, and this is not a small correction if we are 
looking at angular scales smaller than 1/(3. Starting from Eq. (2.7) one can easily compute 
the desired correlation functions. In the case of the angular power spectrum = (a! lm a'l m ) , 
for instance, one obtains 



c\ = ^oj<w {i-\p\e+i+i)^j + v+i)/ 3 ' 



e 3 



3(2* + 1) 



+ W"|£U- (") 



1 Analogous expressions for the CMB Stokes parameters involves spin-weighted spherical harmonics. 
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Thus, for I < 1//3 the leading corrections to the primordial angular power spectrum 
Ci = (a£ m a} m } are small, i.e. 0{(3 2 l 2 ), and moreover they are subject to a further suppression 
since they would give zero in the limit of flat angular power spectrum when Ce ~ C^+i 2 . 
This is the reason why in the literature the aberration has been not considered so far in 
the cosmological parameter probabilistic inference: according to Eq. (2.9) the aberration 
only mildly affects the observed C, and can be therefore safely neglected studying with 
standard MCMC methods the Bayesian credible regions and frequentist confidence levels of 
the cosmological parameters. However as already stressed in the Introduction, this is not 
completely safe for two reasons. First of all, when going at i > the corrections due to 
aberration and Doppler are not under perturbative control and could turn out to be large 
also for the CVs. Second, Eq. (2.9) and (2.8) are correct only in the full sky approximation, 
which does not account for any foreground contamination of the observed maps. In a more 
realistic computation, in fact, one has to consider that indeed only a fraction of the sky is 
experimentally accessible. This can be done by multiplying the original full sky map by an 
opportune window function to mask the background dominated portions of the sky and this 
leads to a larger bias in the Cg%: in fact the small result obtained in (2.9) can be understood 
as a cancellation of a larger effect when summing over the sphere, but such cancellation 
does not hold when we look only at a portion of the sky. More precisely: (1) an effect of 
0(f3l) shows up [6] for asymmetric maps (2) even for symmetric maps the above mentioned 
suppression which acts when Ci ~ Ce+i becomes less effective, because the sum involves a 
larger number of C/s and not only the neighbours. In the present paper we analyze only 
symmetric maps and we leave to further work the analysis with asymmetric maps. 

The procedure of applying a window function leads to a fit of the data based on the 
concept of pseudo-C/. We will review this approach in the next section emphasizing that the 
CMB aberration appreciably modifies the pseudo-C;. This in turn will have a non negligible 
impact on any CMB based cosmological parameter estimation. We conclude this subsection 
mentioning that, because of aberration and Doppler effects, not only the primordial C/ are 
distorted, but also the off-diagonal terms of the covariance matrix (i.e {cL' t ^ m a'l 2m ) with £\ ^ 
£2). In a series of recent work it has been shown that this property of the off-diagonal terms 
can be used to effectively measure the magnitude and direction of our peculiar motion. In 
Ref. [4], for instance, it has been found perturbatively that an experiment such as Planck will 
be able to measure f3 with an accuracy of about 30%, and the direction of this motion with an 
error of about 20 degree, while Ref. [5] shows non-perturbatively that such accuracy would 
be of about 20% for Planck and down to 10% or even 5% for future planned experiments. 
Perturbative corrections to the covariance matrix of the a[ m were also calculated beyond the 
first order approximation. A computation of Wi im a'l 2m ) at second order in f3 can be found 
in Refs. [2, 5], while higher order corrections have be calculated in some cases by [5] and by 
using recurrence formulas in Ref. [7]. 

2.2 Fitting formulas 

Eq. (2.6) can be rewritten in matrix form as follows 

a £m = '^2 K e' e nidi'm ( 2 -10) 



2 Note that such correction goes to zero only if computed with the correct coefficients Eq. (2.7), derived in 
[5] instead of going to 4/3 2 , as previously obtained by [2]. 



- 5 - 



where the kernel Kgi t m is equal to 
and 




s ^sr^i^w- (2 - 12) 



The P™(x) are the associated Legendre polynomials. The K# g m (the so-called aberration 
kernel), which relates the aberrated coefficients a^ m to the primordial ones ai m , are of course 
given at leading order in /3£ <C 1 by expressions such as eq. (2.8), but it is more problematic 
to compute them at large I beyond a perturbative approach. Such problem has been recently 
addressed in [5] proposing an approximate solution to the oscillating integrals of Eq. (2.6). 
In this paper simple fitting formulas were constructed which, approximating the oscillatory 
behavior of the relevant integrals by appropriate Bessel functions, reproduce the exact result 
(verified numerically in a few representative cases) for the elements of the aberration kernel 
with an accuracy of about 0.1 per cent. For example, for the temperature aberration kernel 
an accurate fitting formula is 



Kt-u m ^ Jl{-2c7m) ( 2 - 13 ) 
K e+ Um * M 2 4m) ( 2 - 14 ) 



where J\ is a Bessel function of the first kind. Similar expressions also exist for the polar- 
ization aberration kernels. We refer the reader to [5] for a complete list of fitting formulas 
and a detailed discussion of the derivation and regime of validity. In the next section, we will 
go beyond the approximation schemes discussed in these subsections showing how to include 
exactly aberration and Doppler effects in a MCMC scan of the cosmological parameter space. 



3 Simulated maps including aberration and Doppler 

The ultimate goal of this analysis is to determine the impact of aberration and Doppler on 
a CMB based cosmological parameter estimation. This requires a method to simulate CMB 
maps including aberration and Doppler, in particular when windows functions are considered 
in the analysis to model the CMB maps foreground contamination. Mock data accounting for 
these phenomena can be generated and analyzed by using a straightforward generalization 
of the pseudo-C^ approach. In the following, after introducing the basic formulas to compute 
the pseudo-C^ including aberration and Doppler, we will discuss their implementation in a 
modified version of the HEALPix 3 code [8]. This numerical tool will allow us to generate 
mock data measured by an observer moving with respect to the CMB rest frame. The next 
sections will be then devoted to a MCMC analysis of these simulated data. 

3.1 The boosted pseudo-Q 

To simplify the expression of the aberration kernel (see Eq. (2.11)), we assumed so far that 
the direction of the observer peculiar motion coincides with z, namely the direction of the 
z-axis. However, when dealing with window functions, it is convenient to work using Galactic 

3 http:/ /healpix.jpl. nasa.gov 
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coordinates and assume therefore that the z-axis identifies a direction perpendicular to the 
Galactic plane. We know however that our peculiar motion does not occur perpendicularly to 
the Galactic plane, but it is instead associated with the direction of the CMB dipole, measured 
along the direction identified by the longitude I ~ 249 and the latitude b ~ 48. This implies 
that to express the coefficients computed in Eq. (2.6) consistently with our observed peculiar 
motion, we need to rotate them with a Wigner matrix V mm ,((j),9,j), where <f>,6,j are the 
Euler angles required to align the peculiar velocity, initially along z, with the CMB dipole. 
This Wigner rotation reads as follows 

-e<m'<e 

where the coefficients a' im were calculated in Eq. (2.6). The coefficients ae m are now consis- 
tently expressed in the frame of an observer moving with velocity f3 in the direction of the 
CMB dipole. These equations are correct in a full sky harmonic decomposition. However, 
current CMB experiments observe only a portion of the whole CMB sky. From a theoretical 
point of view, this limitation can be modeled through opportune windows functions which 
set to zero the contributions to the computed quantities coming from regions of the sky back- 
ground dominated (e.g. the Galactic plane). A window function W(fi) acts multiplicatively 
on the original full sky map. Therefore, if we denote by T(n) a full sky temperature map, 
its cut sky counterpart T(n) reads as follows 

r (ft) = w r (fi)f(ft) 

= 2 Yl m m Y im (n)T(n) (3.2) 

l -l<m<l 

where in the second line we expanded the window function in spherical harmonics and denoted 
by W£ m the coefficients of this decomposition. A further expansion of the full sky map 
^( n ) = Ezm^m^m(n) and of the cut sky map T(n) = ^2 tm a em Y £m (n) finally leads to 

l 2 -l2<m 2 <l2 

where the kernel Ti xmx i im2 relating the full sky coefficients a^ 2m2 to the cut sky coefficients 
o< imi admits the following representation in terms of the Wigner 3-j symbols 



* r l\m\l2m2 / j / j w izmz\ l ) 



4-7T 

x (hl2h\(h h h\ {3A) 

\ J \mi —m2 m3 J 

In the present analysis the a,£ m coefficients are the ones given in Eq. (3.1). The coefficients 
d£ m are then used to construct the pseudo-C^, which are denoted by Ci and defined as follows 



Yl i^™i 2 - ( 3 - 5 ) 



-e<m<e 



Combining now Eqs. (3.3) and (3.5), we can finally obtain the desired expression for the 
pseudo-C^ which includes also aberration and Doppler. For the ensemble average of the Ci 
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one finds 



(Ql) — 21 _|_ 1 X/ X/ -^i m ifem2-^7 1 mi^3m 3 (^2m2^l 3 TO 3 ) ■ (3-6) 

—l\<m,\<l\ £2m2 t%ms 

An analytic expression for (C^ 1 ) at first order in f3 has been obtained in Ref. [6]. In the limit 
f3 = and assuming statistical isotropy Eq. (3.6) reduces to 

(C h ) = Y J M h e 2 C £2 (3.7) 
where the mode-coupling matrix M.g x i 2 reads as follows 

^^D-3 +1 )^(r„ 2f „ 3 ) 2 <3 - s) 

and We 3 = 1/(2^3 + 1) ^m 3 l^^"^! 2 ^ s the window function power spectrum. The pseudo-C^ 
formalism is very convenient to construct estimators of the full sky Cg. The best estimate of 
the ensemble average {Cg) is given by C| xp , the value of Cg actually observed or extracted 
from simulations. Then, by inverting the mode-coupling matrix M^^, one can construct 
from Eq. (3.7) an unbiased estimator, denoted by C\, of the full sky Cg 

Ci=E M ll^. (3.9) 

Now, to extract an estimator of the full sky Cg from Eq. (3.6) is instead extremely more 
complicated than in Eq. (3.7), since the boosted Ct involve the full covariance matrix instead 
of only its diagonal part. We will come back to this important observation in section 4 
discussing its role in the cosmological parameter estimation. 

In the usual case of ignoring (3, the goodness of estimators like Eq. (3.9) has been tested 
in various works [12]. In the next sections, when referring to C\, we will use for simplicity an 
estimator obtained replacing in Eq. (3.9) the mode-coupling matrix with the f s k y parameter 
(defined as the portion of the sky covered divided by 4ir). This choice has been shown to 
be fairly robust in the context of weak lensing analyses [13] and it is accurate enough for 
the present investigation, which does not aim at proposing a new estimator to account for 
aberration but rather at determining the error induced by neglecting the peculiar motion of 
the observer when performing analyses based on Eq. (3.9) (with A4^ 2 — > f~ ky )- 

3.2 Numerical implementation 

The equations introduced in the previous sections involve oscillatory angular integrals (see 
Eq. (2.6)) and rotations in harmonic space {e.g. Eqs. (3.3) and (3.1)). To handle these compli- 
cations we will make use of a modified version of the HEALPix code [8] , a numerical package 
to sample functions defined on the sky through a hierarchical, equal area, iso-latitude pix- 
elisation of the sphere. The aim is to simulate CMB data including aberration and Doppler. 
This requires a fiducial cosmological model for which the angular power spectrum Cf is 
assumed to be known. Our choice will be specified at end of this section. The data simu- 
lation is then articulated in five steps which for clarity and completeness we list here in the 
following together with the HEALPix routines used in the generation of the mock data. 



- 8 - 



• synfast: The synfast program generates full sky maps sampling the corresponding har- 
monic coefficients ai m from a Gaussian distribution with zero mean and variance equal 

( f) 

to the Cy . We modified this code including the possibility of applying the transfor- 
mation law given in Eq. (2.5) directly on the pixels. In this way the aberration and 
Doppler effects are included directly in position space rather than in harmonic space. 
As a consequence, one can avoid the evaluation of the complicated oscillatory integral 
in Eq. (2.6). We set the map resolution parameter N s ^ e to the value 2048. 

• anafast: We analyze the temperature map obtained in the previous step with the 
standard HEALPix version of the anafast code, expanding it in spherical harmonics 
and determining the corresponding multipole coefficients. These are the exact (i.e. non 
perturbative) a' em coefficients of Eq. (2.6). 

• alteralm: We also modified the alteralm program in order to rotate these a' im in a new 
reference frame where the peculiar velocity inducing the aberration has the direction 
of the CMB dipole (this transformation is not implemented in the the standard version 
of alteralm). This corresponds to evaluate Eq. (3.1). 

• synfast: Then we performed a second call to the synfast program to generate the 
temperature map T(n) associated with the rotated hi m coefficients. 

• anafast: A final call to the anafast program allows to introduce the desired window 
function and obtain the observed C^ xp . 

This procedure provides the CMB datasets analyzed in the next section. The fiducial model 
used in this paper is similar to the one studied in Ref. [9] . It corresponds to a toy cosmological 
model with five free parameters, namely the baryon density in units of the critical density 
times the square of the present value of the Hubble rate fi^/i 2 , the analogous quantity for the 
dark matter component f^ m /i 2 , the spectral index n s , the amplitude of the primordial power 
spectrum A s and, finally, 100 times the ratio of the sound horizon to the angular diameter 
distance at recombination 0. The fiducial model is characterized by the following parameter 
values: VL^h 2 = 0.022, Sl^h 2 = 0.12, n s = 1, A s = 2.3 x 10~ 9 and has been set to a value 
corresponding to h = 0.7. The other cosmological parameters have been fixed as in the file 
test_params.ini which can be downloaded at the web page [14]. This simplified choice of 
the cosmological parameters allowed us to study the impact of aberration within a model 
for which the Likelihood approximation implemented in the next section has been carefully 
tested in Ref. [9]. 

3.3 Datasets 

Following the procedure introduced in the previous subsection we simulated six datasets, 
full-sky or cut-sky, where an azimuthal cut of 20 degrees around the Galactic plane has been 
imposed. Precisely: two full sky and four cut sky simulations. All simulations will be analyzed 
with Planck-like sensitivity and one of the cut sky simulations also for WMAP-like sensitivity. 
One of the full sky simulations, denoted by SO, assumes that the observer is comoving with the 
CMB ((3 = 0), the second one, instead, includes a peculiar motion (/? = 1.23 x 10~ 3 ) included 
according to Eq. (2.5) and is denoted by SabO. The four cut sky simulations are grouped in 
two pairs corresponding to the same realization of the primordial a£ m , that is, corresponding 
to the same choice of the seed used in the first call to the synfast program (see section 3.2). 
Each pair of "equal seed simulations" consists in one simulation including aberration and 
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Figure 1. Left panel: The red line corresponds to the relative difference between the power spectrum 
estimators of Eq. (3.9) constructed from SI and Sabl plotted as a function of the multipole I. The 
green line is instead the same quantity obtained from the simulations SO and SabO. Right panel: As 
in the left panel but with the cosmic variance in the denominator. 

Doppler, always with j3 = 1.23 x 10 -3 , and one simulation with (3 = 0. The simulations of 
the first pair are denoted by SI and Sabl, being SI the simulation where f3 has been set 
to zero, while the second pair includes the simulations S2 and Sab2, where S2 corresponds 
to the case (3 = 0. In the left panel of Fig.l we show the relative difference between the 
power spectrum estimators of Eq. (3.9) (with M.J^ 2 — > f^} y ) constructed starting from SO 
and SabO (green line) as well as the same quantity obtained from SI and Sabl (red line) 4 . 
We denote this quantity by ACi/Cg. In the comparison between the two full sky simulations 
(SO and SabO) we can appreciate that at large i perturbation theory breaks down and the 
curve oscillates between —0.1 and 0.1. Comparing the two cut sky simulations (SI and 
Sabl), instead, we notice that ACi/Ce is considerably larger than in the full sky case, at 
least up to £ ~ 1000. For larger £, also in this case the curve oscillates between —0.1 and 0.1. 
In the right panel of Fig.l, we finally compare the difference between the power spectrum 
estimators, AC^, with the cosmic variance, which is given by on = ^2/(21 + l)f s k y C\. In 
this figure one can see that for a given I the effect is very large (a few times the cosmic 
variance!). However the difference between the power spectrum estimators oscillates around 
zero without exhibiting clear trends. 

4 For the simulations S2 and Sab2 we only show the parameter reconstruction in section 4. 
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4 MCMC analysis 



Let us assume that a CMB experiment has measured a certain temperature map, say T exp (n), 
with associated angular power spectrum C^ xp . If we wanted to use this map to infer the 
underlying cosmological model, one possibility would be to use directly Eq. (3.9) to derive 
from C| xp the corresponding estimator of the full sky angular power spectrum, C\, and 
implement this quantity in a Likelihood based analysis. This approach would be correct if 
the experiment were comoving with the CMB, but it would instead introduce a "bias" in the 
final parameter reconstruction if the observer were moving with a peculiar velocity different 
from zero, as one can deduce by comparing Eqs. (3.6) and (3.7). The common lore in this 
case is to assume that such a bias is negligible. In the following we will reconsider this 
assumption performing six independent cosmological parameter estimations, one for every 
dataset considered in section 3.3. In these six analyses we will always use the estimator Cj, 
defined in Eq. (3.9) (with MJ^ 2 — > fj^y)-, as an experimental input for our Likelihood based 
reconstruction procedure (described in detail in the next sections). In this way, depending 
from which dataset we are considering, the resulting parameter reconstruction will be biased 
or not (in the sense explained above). Comparing then biased and unbiased reconstructions 
we can estimate the impact of a peculiar motion on the cosmological parameter estimation. 
Summarizing, there are two possible cases: 

• The dataset has been generated including aberration and Doppler. In this case an 
appropriate estimator of the full sky Q could be found solving Eq. (3.7) for the covari- 
ance matrix of the primordial a£ m . Using instead Eq. (3.9) in the Likelihood introduces 
an error in the analysis which consists in interpreting the data with a wrong theo- 
retical assumption, namely assigning a null peculiar velocity to the observer. These 
considerations apply to the datasets SabO, Sabl and Sab2. 

• The dataset has been generated without introducing aberration and Doppler according 
to Eq. (2.5). In this case an analysis based on Eq. (3.9) is perfectly correct. This 
applies to datasets SO, SI and S2. 

Comparing therefore the results obtained for SO and SabO, for SI and Sabl and, finally, 
for S2 and Sab2 will enable us to determine, in one full sky and in two cut sky simulations, 
whether the use of the estimator (3.9), that is, neglecting the peculiar motion of the observer, 
is indeed a completely solid assumption or not. As we will see in this section, this is not true 
a priori in a cut sky analysis. 

4.1 Likelihood 

We will tackle the problem of reconstructing the parameters of our fiducial cosmological 
model from the data simulated in section 3.3 within a Bayesian framework. Bayesian meth- 
ods combined with MCMC scanning techniques allow an effective comparison between ob- 
servations and theoretical predictions. The outputs of this type of analyses are marginal 
posterior probability density functions (pdf) and credible intervals for the underlying model 
parameters (for a review on these subjects, see for instance Ref. [10]). In a Bayesian analysis 
the experimental input are encoded in the Likelihood function, which in our case will be a 
function of the model parameters and of the simulated data. We use here the Likelihood 
approximation introduced in Ref. [9] , according to which the Likelihood C can be written as 



- 11 - 





Figure 2. The top left panel refers to ftbh 2 , the top right panel to fi^m^ 2 , the bottom left to n s 
and, finally, the bottom right to A s . In every panel, the black solid line and the black dot-dashed 
line correspond respectively to the marginal posterior pdf and the mean likelihood derived from the 
dataset SO. The red dashed and dotted lines correspond instead to respectively the posterior and 
mean likelihood obtained from the dataset SabO. In all panels the fiducial model is represented by 
a green dot. In this full sky analysis the statistical outputs obtained from the datasets SO and SabO 
(with and without aberration) exhibit only minor differences. We implemented in the analysis an 
experimental noise corresponding to a Planck-like experiment. 



follows 



- 2 log C ~ Yl 9(Cl/Ct) &P \m^\ ] C\P g(C\,/C fJ 



(4.1) 
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where the function g is explicitly given by 

g(x) = sgn(x - l)y/2(x - ln(x) - 1) , (4.2) 

[M^)] a , is the covariance matrix of the estimator (3.9) evaluated at the fiducial model, C^P 
is the angular power spectrum of the fiducial model and Ci is the predicted angular power 
spectrum (e.g. computed with CAMB). For the fiducial model considered in this analysis, 
this matrix has been already estimated in Ref. [9] and can be downloaded in [14]. In the 
presence of isotropic noise with known angular power spectrum N#, it is straightforward to 
include its contribution to Eq. (4.1). It is in fact enough to sum to all the power spectra 
in Eq. (4.1) a noise contribution Ng. Concerning Ni, we use two noise maps: one which 
reproduces the sensitivity of the WMAP experiment and another for the Planck experiment. 
Regarding the Planck noise we refer to Ref. [9] , implementing in our analysis the same noise 
map used in that analysis. The WMAP noise has been instead taken from [15]. Given a 
signal estimator C\ (see section 3.3), a covariance matrix for it [-W^]^/ an d a noise power 
spectrum Ng, the model parameter reconstruction can be performed with the CosmoMC 
program [11], using these quantities. We restrict our analysis to multipoles t > 30 to avoid 
complication related to the choice of the Likelihoods at small I. Note that, even if those are 
few datapoints, in this way me might actually underestimate the bias we are looking for: in 
fact, as we will see, the power spectrum is more biased due to (3 in a cut-sky at small t rather 
than at large £; however our conclusion that it is not safe to neglect aberration and Doppler 
effect would possibly be only strenghtened in the case in which also the low-£ contribute 
significantly to the bias. Finally, the window function used in this work and applied in the 
last call to the anafast program in the data simulation process (see section 3.2) consists in 
an azimuthal cut of 20 degrees around the Galactic plane. This corresponds to a value of the 
fsky parameter equal to 0.826. 

4.2 Full sky results 

We now show the result of the cosmological parameter reconstruction performed using full sky 
simulated datasets for Planck-like sensitivity. In the four panels of Fig. 2 we plot the marginal 
posterior pdf's and mean likelihoods referring to four parameters of our toy cosmological 
model, namely Ofe/i 2 , Qdmh 2 , n s and A s . For every parameter we show the results obtained 
for both the simulation SO (black lines) and for SabO (red lines) in order to estimate the 
impact of aberration and Doppler on this parameter reconstruction. More specifically, in 
every panel, the black solid line and the black dot-dashed line correspond respectively to 
the marginal posterior pdf and the mean likelihood derived from SO. The red dashed and 
dotted lines, instead, correspond to respectively the posterior and mean likelihood obtained 
from the dataset SabO. From this single map analysis we can conclude that in the case of a 
full sky parameter reconstruction to neglect aberration and Doppler introduces only a minor 
modification to the relevant statistical outputs. Given the null result we only show this for 
Planck-like sensitivity. As we will see, however, this conclusion is not a priori correct in a 
cut sky analysis. 

4.3 Cut sky results 

We now present the results obtained from our cosmological parameter reconstruction per- 
formed with datasets generated including a proper motion of the observer. As in the case of 
a full sky analysis, we report in a single figure the relevant statistical outputs for the param- 
eters Qbh 2 , Qdmh 2 , n s and A s . Fig. 3 refers to an analysis performed using the datasets SI 
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0.985 0.99 0.995 1 1.005 1.01 3.11 3.12 ,„3.13 3.14 

n s log[10 10 A s ] 



Figure 3. As in Fig. 2 but here the datasets analyzed are SI and Sabl. In this figure the means of 
the parameters obtained from SI (no aberration) and from Sabl (with aberration) differ up to almost 
2a (see table 1). The experimental noise corresponds here to a Planck-like experiment. 

and Sabl, while Fig. 4 is associated with the simulations S2 and Sab2, both for a Planck-like 
experiment. Comparing both marginal posterior pdf's and mean likelihoods for the simula- 
tions SI and Sabl we find a more interesting results: the statistical outputs of this analysis 
are biased in a non-negligible way by analyzing the dataset Sabl assuming that the observer 
is comoving with the CMB. We show in Table 1 the difference on best-fit cosmological param- 
eters including or not Doppler and aberration and compare it with the standard deviation for 
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0.99 1 1.01 3.11 3.12 ,„3.13 3.14 

n s log[10 10 A s ] 



Figure 4. As in Fig. 2 but here the datasets analyzed are S2 and Sab2. In this cut sky analysis 
posterior pdf's and mean Likelihoods obtained neglecting aberration (S2) and including aberration 
(Sab2) differ less significantly than in Fig. 3. The experimental noise corresponds here to a Planck-like 
experiment. 

a Planck-like experiment and show that they could deviate up to one or almost two standard 
deviations. In particular the most affected parameters seem to be Q m h 2 and 6. We also find 
a similar effect, though less sizable, for WMAP sensitivity. We show this only for the SI an 
Sabl maps in Fig. 5; the shift in the parameters in this case can go up to one sigma, as shown 
in Table 2. The most affected parameter in this case are given by fl m h 2 and the Amplitude 
of perturbations A s , which are both shifted by one sigma. 
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Figure 5. As in Fig. 3 but with an experimental noise corresponding to the one of a WMAP-likc 
experiment. Also in this case we observe differences between the red (with aberration) and black (no 
aberration) curves which imply a sizable shift in the means of the parameters. In the case of the dark 
matter abundance this shift is close to la. 

Note that this result has a dependence on the specific random seed used to construct 
the map. Actually we have repeated the procedure in a few maps by changing the random 
seed and we found cases where the effect is not so large, as for instance in the dataset Sab2 
(see Fig. 4). Although an analysis performed with a larger number of sky maps, eventually 
using also I < 30, including polarization and considering a more realistic cosmological model, 
is necessary in order to achieve a more robust determination of the impact of aberration on 
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Parameter 


Mean (hi J 


a (SI) 




Mean (babl) 


a (Sabl) 


A (Planck) 


n b h 2 


0.02196 


1.79x10" 


-4 


0.02206 


1.94xl0" 4 


0.60 




0.1194 


1.80x10" 


-3 


0.1168 


1.87xl0~ 3 


1.43 


e 


1.04770 


3.86x10" 


-4 


1.04844 


4.05xl0~ 4 


1.92 


n s 


0.9963 


4.18x10" 


-3 


1.0001 


4.62xl0~ 3 


0.86 


logClO 10 ^) 


3.1281 


4.45x10" 


-3 


3.1230 


4.45x10-3 


1.16 



Table 1. Means and standard deviations for the cosmological model parameters derived by analyzing 
the datasets SI and Sabl using a Planck-like noise. For every parameter we also list in the last column 
the quantity A, namely the difference of the two means (second and fourth column) divided by the 
standard deviation given in the third column. A provides a quantitative estimation of the error 
associated with neglecting aberration effects. 



Parameter 


Mean (SI) 


a (SI) 




Mean (Sabl) 


a (Sabl) 


A (WMAP) 


n b h 2 


0.02146 


3.55x10" 


-4 


0.02166 


3.16xl0~ 4 


0.56 




0.1225 


3.15x10" 


-3 


0.1195 


2.68xl0~ 3 


0.95 


9 


1.04646 


8.88x10" 


-4 


1.04702 


8.08xl0~ 4 


0.63 


n s 


0.9895 


1.05x10" 


-2 


0.9933 


9.35 xl0~ 3 


0.36 


logUO 10 ^) 


3.1386 


6.72x10" 


-3 


3.1319 


5.70xl0-' 3 


0.99 



Table 2. As in table 1 but for the analysis of the SO and Sabl datasets based on a WMAP-like 
noise. 

the cosmological parameter estimation, from the simulations discussed in this paper, we can 
already draw the conclusion that the common lore of neglecting the transformation law (2.5) 
in a CMB based parameter analysis is not in general a solid assumption. 

5 Conclusions 

In this paper we have considered the effect of our peculiar velocity on the CMB maps and 
thus on cosmological parameter estimation. We have assumed a peculiar velocity of f3 = 
1.23 x 10~ 3 c in the direction of the CMB dipole and we have applied a boost transformation 
(eqs.(2.1) and (2.5)) to this frame directly in pixel space on simulated maps, rather than 
on the a£ m 's, which allows us to overcome the challenging problem of evaluating efficiently 
a large number of integrals of highly oscillating Legendre polynomials. For this purpose we 
have simulated a few temperature maps with t < 2100 and with a noise which reproduces 
the Planck satellite resolution as well as maps with WMAP resolution, with £ < 1000, as 
a random realization of a fiducial five-parameters cosmological model, for which a simple 
Likelihood functions has been tested in [9] . For simplicity we have also discarded the first 30 
multipoles, which are usually treated with a different likelihood. Finally we have performed, 
both for full-sky and for cut-sky, an analysis of the maps by running a MCMC on each map 
with or without the aberration plus Doppler effect. 

As a conclusion we find that neglecting aberration would induce some bias on Cosmol- 
gical parameters which could result in a shift of the best-fit values which could go, depending 
on the specific random realization of the map, up to one or almost two standard deviations 
for Planck and up to one standard deviation for WMAP. In particular the most affected 
parameters seem to be O m /i 2 and 9 for Planck. While our analysis is not exhaustive because 
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of many simplifications, we have shown that it is not harmless to ignore aberration in order 
to do cosmological parameter estimation for high resolution experiments such as Planck and 
even for WMAP, and that therefore it would be appropriate before analyzing the experi- 
mental maps to apply a Lorentz trasformation which transforms them back to the CMB rest 
frame. In order to do this it is crucial however to rely on the assumption that our velocity is 
approximately given by the CMB dipole, unless it can be measured to some extent in other 
ways, as proposed in [4, 5]. 
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